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Summary. — The recent measurements of the neutrino reactor angle require a re- 
examination of flavour models based on discrete groups. Indeed, when these models 
deal with the Tri-Bimaximal, the Bimaximal and the Golden Ratio mixing patterns, 
some tensions arise in order to accommodate the reactor angle. In particular, strong 
constraints come from lepton flavour violating processes, like /i — > ey. We present 
the analysis and the main results. 

PACS 14 . 60 . Pq - Neutrino mass and mixing. 
PACS 11.30.Hv - Flavor symmetries. 
PACS 12.60.Jv - Supersymmetry. 



1. Neutrino Data and Predictive Mass Patterns 

Solar and atmospheric neutrino experiments have established the appearance and 
the disappearance of specific flavour neutrinos, that finds the best explanation in the 
oscillation of active neutrinos. Even if the issue of the presence of one or more sterile 
neutrinos must still be clarified, global fits with only three oscillating active neutrinos 
well reproduce the data. Nowadays, the most recent results^) on the oscillation data 
can be summarized in Tab. I. 

More recently, T2K data [1] showed evidences for a non- vanishing reactor angle at the 
3<t level. Subsequently also MINOS [2] and Double Chooz [3] presented their results, in 
agreement with T2K one. In the last months the Daya Bay [4] and RENO [5] experiments 
have released their results on the observation of electron anti-neutrino disappearance, 
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( 1 ) Notice that the best fit values for the reactor angle differ for a factor of 2 in the two fits, due 
to the exclusion (Fogli et al.) or the inclusion (Schwetz et al.) of the data from SBL neutrino 
experiments with a baseline < 100 m. 
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providing at more than 5a and 6c, respectively, the evidence for a non-vanishing reactor 
angle: 

(1) sin 2 9 13 = 0.024 ± 0.005 [Daya Bay] , sin 2 9 13 = 0.029 ± 0.006 [RENO] . 

The average of the results for the reactor angle from the cited experiments, for normal 
(inverted) hierarchy, is given by 

(2) sin 2 6»i3 = 0.022 ± 0.004 (0.023 ± 0.004) . 



Table I. - Fits to neutrino oscillation data. For the fit of Schwetz et all [6], different results 
have been found for the two hierarchies: the IH is shown in the brackets. In both the fits, the 
results correspond to the new reactor fluxes, accounting for the T2K [1] and MINOS [2] data on 
the reactor angle. 
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From the theoretical side, a great effort has been put to construct flavour models that 
are able to describe and explain the experimental results. Before the new data on the 
reactor angle, the attention was focussed on a particular class of mixing patterns, for 
their high predictive power. In the following we will concentrate of the Tri-Bimaximal 
[8, 9] (TB), the Golden Ratio [10, 11, 12, 13, 14] (GR) and the Bimaximal [15, 16, 17, 18] 
(BM) schemes. All these mixing schemes predict a maximal atmospheric angle and a 
vanishing reactor angle, 

(3) sin 2 023 = ^, sin 2 13 = , 

while they differ for the prediction of the solar angle: 

Considering the predicted value of the solar angle for the three mixing schemes, while 
the TB and the GR patterns agree well with the data, the BM one does not at more 
than 5cr. The interest on the BM pattern is mainly due to its relation with the so-called 
Quark-Lepton complementarity [19, 20, 21] (QLC): the QLC consists in a numerical 
relation such that the sum of the experimental values of the lepton solar angle and of 
the Cabibbo angle is roughly 7r/4. From here the idea to revert this expression and write 
6>i2 P w 6»f 2 M - 6 C , where the BM prediction for the solar angle enters [22, 23, 24]. 
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Entering more in details of the these predictive patterns, the unitary matrices corre- 
sponding to the mixing angles listed in eqs. (3) and (4) are the following: 
(5) 

/ 2/^6 1/V3 \ / l/y/2 -l/y/2 \ 

U TB = -l/\/6 l/\/3 -1/V2 , U BM = 1/2 1/2 -l/y/2 , 
V -l/\/6 l/x/3 1/^/2 / V 1/2 1/2 1/72 / 

/ coseg R sm9g R \ 

U GR =\ sm6g R /V2 -cos6g R /V2 l/y/2 . 
V sin^/V2 -cos0f 2 fl /\/2 -l/x/2 / 

In all the three mixing matrices the 13 entry is zero, corresponding to a vanishing reactor 
angle and to an undetermined Dirac CP phase. Moreover, all the entries are pure numbers 
and ensure the independence of the mixing angles from the specific neutrino spectrum: 
this feature is commonly linked to neutrino mass matrices that are form-diagonalizable 
[25] (FD). 

2. Neutrino Flavour Models and the Reactor Angle 

The predictive patterns described in the previous section have been considered as a 
starting point to reproduce the experimental data. To this aim discrete non-Abclian 
flavour symmetries are extremely successful and have been implemented in different ap- 
proaches. In the following we will concentrate on models where the flavour symmetry 
is: global, in order to avoid the presence of a new force, the corresponding gauge bosons 
and their flavour violating effects [26, 27, 28] (a gauged discrete symmetry should be 
considered as a remnant of a gauged continuous symmetry breaking) ; spontaneously bro- 
ken at the high-energy, in order to prevent strong flavour violating effects common in the 
low-energy flavour breaking mechanism [29, 30, 31]; broken by a set of scalar fields, called 
flavons, that transform only under the flavour symmetry, for which a well-defined vacuum 
alignment mechanism can be constructed (this is one of the main advantages with re- 
spect to continuous symmetries [32, 33]). In models that fulfill the previous description, 
the Yukawa Lagrangian is usually written in terms of non-renormalizable operators [34] 
suppressed by suitable powers of the cut-off scale rts A L w A GUT , where A L is the 
scale of lepton number violation and Aqut the GUT scale: 

(6) Cy= A" e i He i + —jy 2AI ' 

Here the Weinberg operators describes the neutrinos, but a completely similar Yukawa 
Lagrangian can be written for the See-Saw mechanisms. When the flavour and the 
electroweak symmetries are broken, the charged lepton and neutrino mass matrices are 
generated. In these models, considering only the lowest dimensional operators, the TB, 
GR and BM patterns could naturally arise as the lepton mixing matrix Upmns- Con- 
sidering also the higher dimensional operators, new contributions correct the LO PMNS 
matrix and are responsible for deviations from the TB, GR and BM predicted mixing 
angles. 

The main ingredient that allows to recover these mixing patterns and their correc- 
tions is the flavour breaking mechanism: the flavons develop vacuum expectation values 
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(VEVs) in specific directions of the flavour space, such that the starting flavour symme- 
try Gf is broken down to two distinct subgroup, G v and Ge in the neutrino and charged 
lcpton sectors, respectively. We will indicate the set of flavons that lead to G v (Ge) as 
$„ G v and Ge represent the low energy symmetries of the neutrino and charged 

lcpton mass matrices: some examples are G v = Z 2 x Z 2 ( 2 ) and Ge = Z n , with n > 3. 

Furthermore, the existence of such symmetry breaking mechanism is usually enforced 
in a supersymmetric context, even if other possibilities have been studied [39, 40]: in the 
following we will consider only supersymmetry flavour models. 

Whether the final PMNS reproduces the experimental data depends on specific fea- 
tures of the models: in the following we identify three major classes that well represent 
the present situation in model building [43, 44, 45]. The GR models can be associate 
with the TB ones for what concerns the results of the present analysis. 

2T. Typical A4 Models for the TB Mixing Pattern. - For this class, we consider for 
definiteness the model in Refs. [39, 41, 42], but the analysis applies to a broader range 
of models based on A4 (see Ref. [43, 44, 45] for details) or on other symmetries (i.e. 
Refs. [46, 47, 48]). The neutrino and charged lepton mass matrices can be written as 

(7) m e = mi a) + Smf 1 , m v = + Sm^ , 

where — diag(y e , y M , y T )vdW, with Vd the VEV of H4 and rj = ($>t)/Af a small 

parameter that breaks down to Ge, and mS> ' is diagonalized by the TB mixing 
matrix. 

In a typical model, the NLO contributions to both the mass matrices correct all 
the entries and are of the same order of magnitude, that we can parametrise with £ = 
($>{;)/ Af rts ($„)/A/, a small parameter that breaks also the subgroups Ge and G v . In 
this case, the mixing angles receive deviations from the initial TB values and we can 
write: 

sin 2 6 23 = \ + Ke(c e 23 ) £ + -L {ne(c^) - V2Ke(c» 3 )) £ 

(8) sin 2 e 12 = l~ \ne{c\ 2 + c? 3 ) £ + ^ Ke(c\ 2 ) £ 



sin 0i3 



3%/2 (c? 2 - cf 3 ) + 2^3 (V2 c? 3 + c^) 



where c^j", complex random number with absolute value of order 1, is the ij entry of 
unitary matrices that diagonalize the charged lepton and neutrino mass matrices at the 
NLO. Accordingly with these expressions, the success rate to reproduce all the three 
mixing angles inside the corresponding 3<r ranges is maximized for £ = 0.07 for both the 
NH and IH. We analyze quantitatively the expressions in eq. (8) and their correlations 
in Fig. 1. The parameters are treated as complex complex numbers with absolute 
values following a Gaussian distribution around 1 with variance 0.5. In the plots we show 
only the NH case. The IH case is similar. 



( ) In some cases, Gf it broken down to G„ = Z2, but an additional accidental Z2 symmetry 
is also present in this sector [35, 36, 37, 38]. 
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Fig. 1. - Typical A4 Models, sin 2 #13 as a function of sin 2 812 (sin 2 #23) is plotted on the left 
(right), following eq. (8). The vertical lines represent the 3a values for sin 2 812 and sin 2 #23, 
following the Fogli et al. [7] fit, in blue, and the Schwetz et al. [6] fit, in red. The horizontal 
lines refers to the 3a values for sin 2 #13 as in eq. (2). 



As we can see, the plots are representing the general behaviour of this class of models: 
sin 2 #13 increases with £, but correspondingly also the deviation of sin 2 #12 from 1/3 does. 
As a result, even for the value of £ that maximizes the success rate, the requirement for 
having a reactor angle inside its 3er error range corresponds to a prediction for the solar 
angle that is no more in good agreement with data. 

2'2. Special A4 Models for the TB Mixing Pattern. - There are some special models 
based on the group A4 [49], in which the LO predictions for the mass matrices are the 
same as in the previous section, but the corrections are not completely generic. In the 
specific case of the model in Ref. [49] , the charged lepton mass matrix still receive generic 
corrections proportional to £ = ($^)/A/, but the neutrino mass matrix is corrected only 
in determined directions: the unitary matrix that digitalize the final neutrino mass matrix 
is given by 



(9) 



U U = U TB V, 



with 



V = 



a 


-e 




where |a| 2 + |£'| 2 = 1, where £' = ($ V )/A/. In this specific model £' > £. The final 
expressions for the neutrino mixing angles after the inclusion of all these corrections are 
given by: 



(10) 
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The success rate to reproduce all the three mixing angles inside their corresponding 
3cr error ranges is maximized by = 0.166(0.171) for the NH (IH). The parameters 
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have been chosen such that £ is a real number in [0.005, 0.06] and cf • are random complex 
numbers with absolute values following a Gaussian distribution around 1 with variance 
0.5. We analyze quantitatively the deviations in eqs. (11) and (12) and their correlations 
in Fig. 2: in the plots on the left (right) column, we show the correlations in eqs. (11) 
and (12) between sin 2 #13 and sin 2 9i 2 or sin 2 #23, respectively: £' is a complex number 
with absolute values equal to 0.166. In the plots we show only the NH case. The IH case 
is similar. For this choice of the parameters, the model can well describe all three angles 
inside the corresponding 3er interval, and its success rate is much larger than that of the 
typical TB models. 




0.15 0.20 025 030 035 0.40 0.45 050 02 03 0.4 03 0.6 0.T 0£ 

Sin 2 *,! SlirSn 
(a) Correlation with sin 2 812 (b) Correlation with sin 2 623 



Fig. 2. - Special A4 Models, sin 2 #13 as a function of sin 2 812 (sin 2 #23) is plotted on the left 
(right), following eqs. (11) and (12). The vertical and the horizontal lines are as in fig. 1. 



2'3. S4 Models for the BM Mixing Pattern. - For the last class, we focus on a rep- 
resentative model based on the $4 discrete group [22]. In this case, mi ' is still the 
diagonal matrix with the charged lepton masses, but is diagonalized by the BM 

mixing matrix. At the higher orders, the neutrino mass matrix preserves the same LO 
flavour structure up to the NNLO level. On the contrary, the changed lepton mass ma- 
trix is corrected at the NLO in all the entries, but not in the 23 and 32 ones. As a result, 
the final neutrino mixing angles at the NLO are given by: 

(13) sme 13 = ^=\ct 2 -ct 3 \^ S m 2 9 12 = 1 --^=ne(cl 2 +cl 3 )^ sin 2 6 23 = \ . 

To properly correct the BM value of the solar angle to agree with the data, £ is expected 
to be O(Ac). Studying the success rate to have all the three mixing angles inside the 
corresponding 3cr ranges, we find that it is maximized for both the NH and IH when 
£ = 0.163. We analyze quantitatively the expressions in eq. (13) and their predictions 
in Fig. 3, where 012,13 have been taken as random complex numbers with absolute value 
following a Gaussian distribution around 1 with variance 0.5, while £ = 0.185(0.194). 
A value close cos 6 close to —1 is favoured in order to maximize the success rate 
[43, 44, 45]. In Fig. 3, only for the NH case is shown. The IH case is similar. 

Considering the results for the success rates of all the three classes of models, these S4 
models arc strongly disfavoured with respect to the special S4 ones, while are comparable 
with respect to the typical A4 models. 
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Fig. 3. - S4, Models, sin 2 #13 as a function of sin 2 #12 is plotted, following eq. (13). 
3. Conclusions 

Discrete symmetries can well accommodate the neutrino mixing pattern, especially 
considering an approach in which the PMNS matrix is given in first approximation by the 
TB, the GR or the BM patterns. With the new results on the reactor angle, however, it 
appears suspicious that one of these mixing schemes could be a fundamental structure of 
nature, while it is getting stronger the feeling that they are simply numerical accidents. 
Indeed, the type and the size of the corrections necessary to bring these mixing patterns 
in agreement with the data put sever doubts on their naturalness. 

Furthermore, as discussed in a series of papers [50, 51, 52, 53, 54, 55] and updated in 
Ref. [43, 44, 45], the analysis of lepton flavour violating transitions is fundamental to test 
flavour models. In particular, with the new results on the reactor angle and the large 
size of the NLO corrections, the bounds on the supersymmetric parameters space coming 
from the ti — » £jj decays are strong, even for small tan/3, and if light supersymmetric 
particles are found then these models are disfavoured. Moreover, it appears impossible to 
satisfy the MEG bound and, at the same time, to reproduce the muon g — 2 discrepancy. 

Even though the huge effort of these years in constructing flavour models to describe 
masses and mixings for the neutrinos, and more in general for all the fermions, it is dis- 
couraging that no illuminating strategy arise form this scenario. On the other hand, this 
is partially related to the large uncertainties still present in the flavour sector. The hope 
is that with a better determination of the lepton mixing angles and with the knowledge 
of the CP phases, the neutrino mass scale, the type of the neutrino nature and spectrum, 
it will be finally possible to shed light on the origin of the fermion masses and mixings. 

* * * 
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